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Abstract. We investigate a new representation of general operators by means 
of sums of shifted Gabor multipliers. These representations arise by studying 
the matrix of an operator with respect to a Gabor frame. Each shifted Gabor 
multiplier corresponds to a side-diagonal of this matrix. This representation 
is especially useful for operators whose associated matrix possesses some off- 
diagonal decay. In this case one can completely characterize the symbol class 
of the operator by the size of the symbols of the Gabor multipliers. As an 
application we derive approximation theorems for pseudodifferential operators 
in the Sjostrand class. 



1. Introduction 

One of the approaches to understand a given operator is to decompose (= "an- 
alyze") the operator into simpler operators and then to find approximation the- 
orems. The meaning of "simple" and "approximation" varies with one's point of 
view and with the application at hand. In this investigation we take the point of 
view of time-frequency analysis and consider an operator as simple if it aligns well 
with Gabor frames. Technically, the simple operators are Gabor multipliers. Given 
a point z = (x,£) G M? d in the time-frequency plane, we denote the associated 
time-frequency shift acting on a function / by 

n(z)f(t) = e 2 ^ t f(t-x), x,^teR d . 

Now fix a non-zero function g G L 2 (IR d ) (a so-called "window" function) and a 
lattice A = AZ* 2d (with an invertible 2c? x 2d-matrix). A Gabor multiplier with the 
symbol a = (g^aga is defined informally as 

(1) M^/ = £a A </,*(A)<7>7r(A)0. 

AeA 

Gabor multipliers have attracted attention both in engineering and in mathemat- 
ics [2,5,9,10,23], because they provide an easy and computationally attractive 
model of time-frequency masking. In ([I]) the Gabor coefficients (f,7r(X)g) are a 
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measure of the time-frequency content at A = (Ai, A2) G M? d , i.e., the amplitude of 
the frequency A2 at time Ai. The multiplication of the Gabor coefficients with the 
symbol or mask a amounts to an enhancement or damping of certain time-frequency 
regions. In this sense, the symbol a can be compared to the transfer function of a 
discrete time-invariant system usually given by a convolution operator f —>■ f * a. 

Whereas the basic properties (boundedness, mapping properties, Schatten class 
properties) of Gabor multipliers are well understood [10], it is still open how Gabor 
multipliers fit into the general picture of (pseudodifferential) operators. How big is 
the class of Gabor multipliers? How well can a given operator be approximated by 
a Gabor multiplier or by a sum of Gabor multipliers? Which properties influence 
and determine the accuracy of an approximation by Gabor multipliers? 

To address these questions, we will prove results of the following type. 

(1) We show that every operator A from iS(lR d ) to iS'(IR d ) can be represented as 
a sum of shifted Gabor multipliers 



for suitable windows g and symbols sequences a„. This decomposition is a new 
representation of operators comparable to the Schwartz kernel theorem (which we 
will indeed use to deduce (J2])). 

(2) We relate the properties of the symbol of a pseudodifferential operator and 
the multiplier sequences sl u . One of the main results is a characterization of the 
generalized Sjostrand classes by the size of the multipliers a u . 

(3) We estimate the error 



in the operator norm on L 2 (IR d ) and on other function spaces. In this way we 
develop an approximation theory for operators. One can draw an analogy to the 
classical theorems of Jackson for the approximation of a continuous function by 
trigonometric polynomials, e.g., [3]. Instead of continuous functions we consider 
bounded operators, and instead of trigonometric polynomials we approximate with 
finite sums of Gabor multipliers. Our main result is in the spirit of Jackson's theo- 
rem: the rate of approximation by sums of Gabor multipliers is directly correlated 
to the smoothness of the symbol. 

(4) Finally we indicate how Gabor multipliers are in the modeling time-varying 
environments in wireless communication. 

A related investigation of Gabor multipliers was carried out by Dorfler and Tor- 
resani [5]. They approximate Hilbert-Schmidt operators by Gabor multipliers and 
sums of Gabor multipliers with different windows and perform a window optimiza- 
tion. Though similar in spirit, the results and methods are completely disjoint. 
The Hilbert-Schmidt norm of an operator is accessible to explicit calculations, and 
questions of best approximation and error estimates can be treated by using or- 
thogonal projections. By contrast, there is no formula for the operator norm of an 
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operator, and we have to use different ideas. Our results are based on the almost di- 
agonalization of pseudodifferential operators with respect to Gabor frames [13] and 
the analogy between matrices and operators. In this analogy a Gabor multiplier 
corresponds to a diagonal matrix, and the representation (jSJ) corresponds to the 
decomposition of a matrix into the sum of its side-diagonals. The approximation 
in the operator norm is also motivated by the needs of wireless communications. 
There the arising operators are perturbations of the identity operator and of con- 
volution operators and they can never be Hilbert-Schmidt operators. 

The paper is organized as follows: In Section 2 we recall a minimum of definitions 
and results from time-frequency analysis, in Section 3 we introduce the main tool for 
the approximation theory by Gabor multipliers, namely the almost diagonalization 
of pseudodifferential operators in the Sjostrand class. Section 4 presents the formal 
definition and boundedness properties of Gabor multipliers. In Section 5 we derive 
several versions of the operator representation (j2j), and in Section 6 we investigate 
the approximation of operators by Gabor multipliers. We conclude with a brief 
discussion of Gabor multipliers in wireless communications in Section 7. 

2. Some Time-Frequency Analysis 

Let z = (x, £ R d x be a point in the time-frequency plane, t G M. d . The 
time-frequency shift ir(z) is the operator 

n{z)f{t) = e 2 ^ t f{t-x) 

acting on L 2 (M. d ), S'(M. d ) and many other spaces. In the following we use [11] as a 
general reference to time-frequency analysis. 

Associated to this set of operators are a signal transform, function spaces, and 
structured frames; namely the short-time Fourier transform, the modulation spaces, 
and Gabor frames. 

Fix a non-zero window function g, g G S, say. The short-time Fourier transform 
of / with respect to the "window" g is defined as 

VJ(z) = (f, n(z)g) = [ f^W^e 2 ^ 4 dt = (f ■ g(. - x))^) , 

whenever the duality (•,•), the integral, or the Fourier transform are defined. 

Further, we may associate a class of function spaces to the time-frequency shifts 
and the short-time Fourier transform. Again for fixed non-zero window function g, 
g G <S, say, we define the modulation spaces M^f for 1 < p, q < oo and m a weight 
function as follows: 

/ G M%f(R d ) & V g f G 

with norm = II^/IIl^- Here L p ^ is the usual mixed-norm space. These 

spaces are well-defined, there is an extensive theory about modulation spaces. The 
reader should consult [11, Chps. 11-13] and the references therein, an important 
source for the history is [6]. 

Another notion involving the set of time-frequency shifts tc(z) is the notion of 
Gabor frames. Fix a lattice A C ]R 2a! , i.e., A = AZ 2d for an invertible 2d x 2d 
matrix A. The collection of time-frequency shifts Q(g,A) = {7r(X)g : A G A} for 
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some non-zero g G L 2 (IR d ) is called a Gabor system. The set Q(g,A) is a Gabor 
frame, if there exist constants A, B > such that 

(3) A\\f\\t < J2\(f^(\)g}\ 2 < B\\f\\l V/ G L 2 (M. d ) . 

AeA 

If A = B = 1, Q(g, A) is called a Parseval frame and (J3J) implies the expansion 

(4) / = X;</,ir(A)«/MA)0 V/GL 2 (M d ) 

AeA 

with unconditional convergence in L 2 (R d ). In contrast to an orthonormal expansion 
the coefficients in this expansion need not be unique, the choice (/, it(X)g) is a 
distinguished, explicit, and convenient choice. The existence problem for Gabor 
frames is almost completely settled [1]. In what follows, we take the existence of 
Gabor frames and Parseval Gabor frames for granted. 

Characterizations of modulation spaces by Gabor frames. In the fol- 
lowing v always denotes a submultiplicative, even weight on R 2d , i.e., v(x + y) < 
v(x)v(y). A weight function m on M? d is called w-moderate, if m(x+y) < Cv(x)m(y) 
for all x, y G M. 2d . 

If Q(g, A) is a Gabor frame and g G M^(R d ), then magnitude of the frame coeffi- 
cients (/, 77(A)(7), A G A, characterizes the membership of a function in a particular 
modulation space [8,11]. More precisely, for every p G [l,oo] and v-moderate 
weight m, 

f G M%f(R d ) & ( I (/. < X )9) \ p m(\y) UP < 00 , 
AeA 

and the latter expression is an equivalent norm on M^f. An analogous characteri- 
zation holds for rectangular lattices and the mixed modulation spaces [8]. For 
p 7^ q and a general lattice A, we have / G M™, if and only if J2\eA I (/' 7r (^)fi f ) Ixa+q G 
L^(R 2d ) [7]. ^ ' 

In particular, we have the following simple characterization of the Schwartz class 
by means of Gabor frames. Assume that g G and that Q(g, A) is a frame for 

L 2 (R d ). A function / G L 2 (R d ) belongs to S(R d ), if and only if (/, ir(\)g),\ G A 
decays rapidly, i.e., 

(5) |</,7r(A)flr>| = 0((1 + |A|)- S ) for all s > . 

Pseudodifferential Operators. Given a function (or distribution) a on M 2d , 
the corresponding pseudodifferential operator with symbol a is defined informally 
by the integral 

(6) a{x,D)f{x)= f a(x,0f(0e 2mx< d£. 

Jm. d 

Again the definition of pseudodifferential operators does not necessarily require 
that the integral is defined, see [12, 16] for a rigorous definition. 
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Symbol Classes. We consider certain modulation spaces as symbol classes. 
They are generalizations of the Sjostrand class [18] and may be understood as non- 
smooth extensions of the Hormander class Sq . Fix a non-zero test function $ 
on M. 2d , e.g., the Gaussian $(V) = e~ nz ' z and a submultiplicative weight v on M. 2d . 
Then a symbol a belongs to M^'^M 2 ^), if 

(7) / sup \V$a(z,()\ v(()d(= \\a\\ M ~,i < oo . 

Jm. 2d z m 2d 

Likewise, a symbol a belongs to M^'°°(K 2d ), if 

(8) sup sup | V$a(z, ()\v(() = \\<j\\ms°-°° < 00 • 



The generalized Sjostrand classes are thus certain modulation spaces on M. 2d . 
Further generalizations use the class of solid convolution algebras on J? d to param- 
etrize the extensions of Sq [14]. 

3. Almost Diagonalization 

We will study the matrix of a pseudo differential operator a(x, D) with respect 
to a Gabor frame. 

Precisely, let Q(g,A) be a Gabor system and a a symbol with corresponding 
pseudodifferential operator a(x, D). Then we consider the matrix M = M(a) with 
entries 

(9) M(a) v = (a{x, D)(n(fj)g),ir{X)g) X, fi E A . 

The main results of [13, 14] establish a precise link between pseudodifferential 
operators in generalized Sjostrand classes and the associated matrix. In the follow- 
ing j denotes the rotation j(zi, z 2 ) = (^2, —Zi). (It is needed to use the standard 
modulation operators instead of some symplectic modulations.) 

Theorem 1 ( [13]). Assume that g G M^(IR d ) and that G(g,A) is a Gabor frame. 
Then a G M^^Li if and only if there exists h G ^J(A) such that 

(10) \M(a) x J = \ (o-(x,D)n(^g,n(X)g}\ < h(X-fj) 

for all A, G A. Furthermore, inf \\h\\ti with the infimum taken over all h satisfying 
([10]) is an equivalent norm on . 

REMARK: If G(g,A) is a merely a Gabor system with g G M] .-i(R d ) (but not a 
frame), then the almost diagonalization ( flOl) still holds. However, in general, the 
converse is not true in general, because the matrix M(o~) does not fully describe 
the operator a(x, D). 

A similar theorem holds for the class ill? ' 00 . 



V 



Theorem 2 ( [14]). Assume that g G M^,(R d ) and that Q(g,A) is a Gabor frame. 
Then a G M™'^ if and only if there exists 

(11) \M{a) x J = I (<t(x, D)7r(fj,)g, ir(X)g) \ < Cv(X-ft)- 1 

for all X, // G A. 
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Thus in the case of M^ ' 00 one obtains genuine off-diagonal decay of the associated 
matrix M(a). The standard weights to describe the decay condition are either 
polynomial weights v(z) = (z) s = (1 + |z| 2 ) s / 2 or (sub)exponential weights v(z) = 
e a|2|b for a > and < b < 1. 

Other types of decay conditions are studied in [14]. 

As a special case we mention the Hormander Class Sq consisting of all C°°- 
functions on R 2d with bounded derivatives. This symbol class is related to modula- 
tion spaces as follows [14,22]: Sq >0 = f] s>0 Mj^f°. As a consequence of Theorem [2] 
we obtain a characterization of Sq Q . 

Corollary 3. g G S, g ^ 0, Q(g,A) Gabor frame. Then a G Sq if and only if for 
every s > there is a constant C s such that 

| (a(x, D)ir(ji)g, n(X)g)\ < C s (1 + | A - fi\)- s for all X, p. G A . 

For completeness we mention that operators with symbols in one of the gener- 
alized Sjostrand classes M^ ' 1 and M^ ' 00 are closed under composition and that 
they are bounded on many modulation spaces [13, 14, 19]. 

4. Gabor Multipliers 

Next we introduce a special class of operators, so-called Gabor multipliers. These 
operators are particularly simple and, in some sense, correspond to diagonal ma- 
trices. 

Definition 1. Let Q(g,A) be a Gabor system. Given "symbol sequence" a = 
(oa)a6A) the Gabor multiplier M a is defined to be the operator 

M a f = a M> 1T (p)g}'rr(p)g . 

Clearly, the definition also depends on the Gabor system Q(g,A). To indicate 
the dependence of the Gabor multiplier of all parameters, we would have to write 
M|' A . Since our results holds generically for all Gabor systems in a certain class, we 
prefer to keep the notation simple and omit the reference to Q(g, A). Likewise one 
could use two windows g and 7 and consider Gabor multipliers of the form M a f = 
S M eA a A*(/' 71 '(/-0 d) ^ (t l )l- Such generalizations cause only notational complications, 
but do not change any of the results or arguments. 

Gabor multipliers have been studied in detail by Feichtinger and Nowak [10]. 
Gabor multipliers can also be interpreted as time-frequency localization operators 
with distributional symbols, and boundedness results follow from the theory of 
localization operators [4]. 

As a typical boundedness result we state the following one. 

Lemma 4. Assume that g G Ml (IR d ) . If a G £°°(A) ; then M a is bounded on every 
modulation space for 1 < p, q < 00 and every v-moderate weight m with a 
uniform bound for the operator norm 

(12) ||M B || <c A |MI M 1 ll a l|c» • 
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Proof. The decisive property for the proof is the boundedness of the coefficient 

operator / — > ( (/, Tc(X)g) ) from M^fiW 1 ) to the sequence space -^ 9 (A) and 
V / AeA 

the synthesis operator c — > J2 X eA c ^ 7T W9- See [H> Thms. 12.2.1-4] or [8] (for 
rectangular lattices). Using these estimates, we find for p = q that 



\MJWug* 



^2a^(f,Tv(n)g)7r(fi)g\\ M P, P 



fieA 



< 



C\\g\\ Mi (J2K(f^^)9)\ P m^T) 



Vp 



M6A 



(13) 



< C\\g\\ Mi ||a|| 00 (^|(/,7r(^)rm(^) P ) 1/P 

< (C|b|| M i) 2 ||a|UI|/IU^- 



For p ^ q and non-separable lattices the proof is identical, once the sequence space 
is defined correctly. See [7]. ■ 

5. Representations of Operators with Gabor Multipliers 

After these preparations let us now explain how and why Gabor multipliers 
arise in the theory of pseudodifferential operators. Assume that Q(g,A) is a Par- 
seval frame, so that (j3J) holds. We expand both / and Af = a(x,D)f with 
respect to the frame Q(g,A). Then / = ^)„ e ^(/, 7r (fi)g)7i(fi)g and a(x,D)f = 
^ AgA (v4/, n(X)g)7r(X)g. Substituting the expansion of / into the coefficients (Af, 7r(X)g), 
we obtain that 

Af = ]T(Af,7r(A)s>7r(A)s 

AeA 

AeA /^eA 

= Yl (^2(^(^)9, tt(^ + v)g) (f, Tr(n)g)^ n(fi + v)g 

u&A fj,eA 

(14) = < v ) (j2(Mri9, ttOi + u)g) e 2 ™^ (/, nQi)g) n^)g) 

ueA fj,eA 



In the transition to the last line we have used the commutation rule for time- 
frequency shifts 7r(/x + v) = e 27rWl M2 7r(^)7r(/x). 

This calculation is known, of course, to every student of linear algebra; it shows 
how to express a linear operator by its matrix with respect to a "basis" , which in 
our case we take to be a Gabor frame. 

Now define the sequence a u by 



(15) 



a„(//) = (An(n)g, tt(/x + v)g) e 
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then the expression in (j!4p in parenthesis is just a Gabor multiplier with symbol 
a u , and we may rewrite ffT4l) as 

(16) a(x,D) = J2Av)M a „. 

i/eA 

At least informally, the above identity shows that every pseudodifferential operator 
is a sum of shifted Gabor multipliers. In this representation each Gabor multiplier 
corresponds to a side-diagonal of the associated matrix M(a). 

To make this argument precise, we offer several versions. We first show that every 
reasonable operator can be represented as a sum of shifted Gabor multipliers. 

Proposition 5. Let Q(g,A) be a Parseval frame and g G iS(M d ) \ {0}. 

(i) If the sequences sl u satisfy the growth conditions 

M/i)l < C(l + \n\ + \v\) N for A, f^eA 

then ^2 ueA 7i(i/)M av defines a continuous operator from S(R d ) to S'(M. d ) and the 
series converges in the weak operator topology. 

(ii) Conversely, assume that A is a continuous operator from S(M. d ) to <S'(R d ) 
(with the weak* -topology) . Then A possesses a representation (|T6l) as a sum of 
shifted Gabor multipliers with symbols a„ satisfying the growth estimate 

\aM\<C(l + \fx\ + \u\) N 

for constants C, N > depending only on g and A. 

Proof, (i) Fix v G A and consider the Gabor multiplier M &v f = 
= ^2 eA a. u (fi)(f, Tr(n)g)ir(/j,)g. Since by hypothesis a u grows polynomially and 
(/, ft((J>)g) decays rapidly by (jSJ), the coefficients a I/ (//)(/, Tr{n)g) also decay rapidly, 
whence the sum converges in S(R d ) and M a J G S(R d ). Consequently the partial 
sums Y1i\v\<m M av f are in 5(lR d ), and we only need to verify that the partial sums 
converge in the weak* sense to an element mS'(R d ). The convergence follows from 

\u\<M \v\<M 

< ct E^+iH + inra+iHr" 

\u\<M 

for arbitrary I > 0. Hence the sum converges, or equivalently, the infinite sum 
J^veA-Ma^f converges in the weak-* sense. 

(ii) Let A be continuous from 5(R a! ) to S'(M. d ). As suggested by informal deriva- 
tion (TH1) . we choose the symbols a u to be a y (//) = (A7r(fi)g, 7r(/i + v)g) e 27rWl '^ 2 . 

By the Schwartz kernel theorem there exists a tempered distribution k G iS'(!R 2d ), 
such that (Af, /i) K d = (k, /i<g)/} R 2d for /, h G <S(M d ). Since the STFT of a tempered 
distribution grows at most polynomially, there exists an N > 0, such that 

\(k,ir(\)g®n(n)g)\ = |V ff ® s fc(Ai, m; A 2 , -/j, 2 )\ < C(l + |A| + \fi\) N VA,//G A. 
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We deduce the claimed growth estimate for the symbols as follows: 

|a„(/i)| = \(Air(iJ,)g,n(ji + u)g) Rd \ 

= \(k,7r(ij + v)g (g)Tc(fi)g) R 2 d \ 

< C(l + \»\ + + u\f < C"(l + + \u\) N . 

By Step (i) the sum J2ueA 7r ( z/ )^a I/ is a continuous operator from 5(lR d ) to 
5'(R d ), and by the choice of the multipliers a u this operator must coincide with 
the given operator A. ■ 

Our next version of ( fl6l) deals with the representation of pseudodifferential op- 
erators. 

Lemma 6. Assume that Q(g,A) is a Parseval frame and g G M^(M. d ). 

(i) If the sequence of symbols a„ satisfies the condition 

22 \\a-v\\oov{v) < oo , 

i/GA 

then the sum of shifted Gabor multipliers ({16]) converges in the operator norm on 
for every 1 < p, q < oo and every v-moderate weight m. 

(ii) If a G M^j.j (IR 2d ) ; then the series ( fl6l) converges in the operator norm of 
to cr(x,D), and X^eA H^IU^) ^ ^INI a/ 00 * 1 , ■ 

Proof, (i) By ( fl2l) the operator norm of M &v on is bounded by CHa^H^. 

Further, the operator norm of a time-frequency shift 7r(z/) on M^f is bounded by 
v(i/). Thus, summing over u, we obtain 

\W(x,D)\\ M p,*^ M P l i < ^\\Tr(v)M av \\ M p,*_. M ™ 

ueA 

ueA 

< 22 ll a ^ll°° v i v ) < 00 • 

Thus the series of operators (fl6l) converges absolutely in the operator norm on 
MP'i. 

m 

(ii) First note that by Theorem [1] there exists a h G ^(A) such that \\h\\ti < 
C\\a\\ M oo,i and 

= |(a(x,D)7r(^,7r(^ + %) e 27 ™| < %) 

By Step (i), ffl6|) converges absolutely in the operator norm on M^' 9 , and by the 
choice of symbols the limit coincides with a(x, D). The norm estimate follows from 
EkeA \\ a v\\oov(v) < \\h\\ e i < C\\a\\ M ^i_ * 

voj 

REMARK: By adjusting the conditions in the above proof, many more versions of 
({16]) can be shown to be meaningful. 
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By combining the almost diagonalization of Theorems [I] and [21 we now charac- 
terize operators with symbols in the generalized Sjostrand class by means of sums 
of Gabor multipliers. 

Proposition 7. Fix a submultiplicative weight v and a non-zero g G M^(M. 2d ) such 
that Q(g,A) is a Pars eval frame. 

(i) A symbol a belongs to M^j_ 1 (R M ), if and only if there exist sequences a„ G 
£°°(A), such that 

(17) a(x,D) = J2^)M a „ 

and 



Furthermore, the sequence of multipliers a„ can be chosen such that^2 ueA \\b. u \\ 00 v(i') < 
C||<x|| M oo,i < C'E^eA hv\\ooV{v). 

(ii) Assume in addition that t> _1 * v~ l < Cv~ x (v is subconvolutive) . Then 
a G M™'™! (M. 2d ), if and only ifa(x,D) possesses a representation (fT7|) with 

||a„||oo < Cv(v)~ l . 

Again the multipliers a„ can be chosen, such that sup u£A WsluWooV^) is an equivalent 

Moo,oo 

VOJ 1 

Proof. The sufficiency was shown in Lemma [61 If o G then both f[T8|) and 

f[TT|) are satisfied. 

Assume conversely that an operator A : iS(lR d ) — > S'(M d ) is given as a sum of 
Gabor multipliers (fTT|) with symbols satisfying f[T51) . Then by the Schwartz kernel 
theorem A possesses a symbol a G iS'(IR 2d ) and A = o~(x, D). 

To show that a is in M^° ,1 (]R M ), we estimate the size of entries M(a)\^ with 
respect to a reference frame Q(<f, A) and then apply the characterization of Theo- 
rem [TJ 

Fix a Gabor frame Q(if,A) with ip G 5(M d ) H M^(M. d ). Then the matrix entries 
M(a)\fj, = (a(x, D)(7r([j,)<p),7r(\)(p) are well-defined. Since ^2 ue \M a „ converges 
weakly, we may interchange the brackets (•, •) with the summation over v and 
obtain 

M(a) Xfl = (o~(x, D) 7r(/i)(/?, 7r(A)(/?) 

(19) = ^Tm^m^v^av) 

= "£2^a-v(K)(7r({i)ip,7r(K)g} (Tr(y)ir(K)g,ir(X)<p) . 
^eA kgA 

Since g G Af* and G 5(M d ) n M*(IR d ), the sequence 

h(\) = \{g,n(\)<p)\ 
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belongs to ^J(A) and \\h\\gi < C\\g\\M^ (with a constant depending only on (p and 
A) [11, Prop. 12.1.11]. Set a(u) = Ha^Hoo, h*(X) = h(—X), take absolute values in 
f fl9|) . and use the decay of h. Then we obtain that 

\M(a) XtJ \ < y~] ^ ll^lloo h{ji - k) h(X - v - k) 
= ^ a(u) ^ h*(K)h(X - v - {k + n)) 

= ^2a{u)(h**h)(X-fjL-u) 

ueA 

= (a * h* * h)(X — //) . 

Since a G £l by assumption (|T7|) and h E l\ because g G M*, the matrix M(cr) is 
dominated entrywise by the sequence a * h* * h G i\. 

Theorem [1] applies and we conclude that a G M^_ 1 (R 2d ), as claimed. Further- 
more, 

IMIm 00 ' 1 < C|| q; * /i* * h\\ii 
< C||q||^||/»||^ 

The proof for a G M^j^ is similar. The only modification occurs in the last 
part, where we have to use the subconvolutivity of v and the convolution relation 
C * C C C ■ 

6. Approximation Theorems 

Since every operator A can be represented as an infinite sum of shifted Gabor 
multipliers, it is natural to truncate the infinite series (IT7|) and approximate A by 
a finite sum of Gabor multipliers. Compared to a general operator, a finite sum of 
Gabor multipliers is easy to understand and easy to treat computationally. 

For the formulation of the approximation theorems we introduce the error 

(20) E N (a) := \\a(x, D) - J] n{v)M a „ \\ m ™-*m™ ■ 

\v\<N 

We take some liberty in the interpretation of the operator norm involved. As we 
have seen, the spaces we can chose depend mostly on the quality of the window. 
Here are some precise approximation theorems. 

Theorem 8. Assume thatQ(g,A) is a (Parseval) Gabor frame for L 2 (M. d ) . 

(i) If g G M 1 (M d ) and a G M 00 ' 1 , then Ejy(a) — ► (in the operator norm on 
for alll<p,q<oo). 

(ii) Ifge Ml and a G M^Z, then 

En(o~) < C IHIm 00 ' 00 / v { u ) 1 ■ 

\v\>N 
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(in) If g G Ml and a G M^i, then 

En(ct) < C Hcrllj^oo,! sup v (z^) _1 . 



v °i 1 \u\>N 

(iv) If g G M^ s and a G M^' s °° and s > 2d, then 

? j N (o ) \ O || O || M 



Efj(a) < C Hull m°°'°° -/V 2d 



'(0 s 

(v) If g G ML S and cr G M^ 1 , then 

E N (a) < C\\a\\ M oo,i N~ s . 

Proof, (i) By Proposition [7| we represent a(x, D) as a sum of Gabor multipliers 

a(x,D) = Y^<v)M av , 
ueA 

such that ^ vgA ||a-j/||oo < °°- Since the sum converges in the operator norm on 
M p,g , the difference between a(x, D) and its approximation becomes 

a(x,D) - 5>(»M a „ = £ n(v)M a „ . 

i/£A \u\>N 

Taking operator norms, we obtain the following estimate for the error En(o~): 
E N (cr) = \\a(x, D) — ix{y)M Av \\ M p, q ^ MP , q 

= II ^2 n ( l/ )M aL J\ M p, q ^ M p, q 

\u\>N 

— H^ 3 " l|MP.9-»MP.9 

M>AT 

(21) < ^ KHoo. 

M>AT 

Since a G M 00 ' 1 and J2ueA ll a ^ ll°o < oo, (121]) implies that En (a) — > 0. 

(ii) By Proposition[7|the assumption a G M^.^ implies that Halloo < C||a"|| M °°.°° 

Consequently, by (I2TI) the approximation error is at most 

E N {cr) < V" lla^Hoo < C||cr|| M ococ V" v^)' 1 . 

\u\>N \v\>N 

If = (1 + |z|) s , then ^ij,i >Ar (l + M)~ s < CN 2d ~ s , whence assertion (iv). 

(iii) Likewise, if a G M^.^, then by Theorem [T] 

(u) < C||cr|| M oo,i < oo 

veA 

and thus 

E N {cr) < ^2 H a «H°o ^ SU P v ( v Y l ^2 

\u\>N H>N H>N 
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If v(z) = (1 + \z\) s , then supi^jy v{y) 1 < (1 + N) s , whence (v). 



REMARKS: 1. The error estimates of Theorem [8] yields an estimate for the size 
of the cut-off parameter N and thus of the number of shifted Gabor multipliers 
required for a good approximation. For instance, if a G M,^' s (R 2d ) and a tolerance 
e > are given, then the error estimate 



E N (a) < C |M|. f0 c,i AT 



implies that 



e 



JV > ' " 10 



It is possible to say more about the dependency of the constant C on g and A by 
precise bookkeeping in the results about almost diagonalization of pseudodifferen- 
tial operators in [13]. 

2. One may wonder whether a converse of Theorem [S] holds and whether the qual- 
ity of approximation as expressed by Theorem [H] characterizes the symbol classes. 
This guess is false, because the operator norm \\a(x, D)\\m™-+mp>* an d HcIIm 00 - 1 are 
not equivalent. Indeed, following an idea of Klotz [17] about approximation alge- 
bras of matrices, one can define a new symbol class directly by the approximation 
properties of Gabor multipliers as follows: We say that a G .4 s (lR 2a! ), if 

inf \\a(x, D) — V M a J L ^ L 2 < CN~ S . 

a„e^°°(A) 

v&A,\v\>N 

This class of symbols contains M^' s °°. By comparison with the corresponding 
matrix algebras one can verify that M^' s °° is strictly smaller than ^4 s (IR 2<:i ). 

7. Application 

Finally we mention the problem that has motivated the approximation theorems 
in the previous section. This problem concerns the transmission and decoding of 
digital information by a variant of (orthogonal) frequency division multiplexing 
(OFDM). 

Here is a very coarse description of this procedure. 

1. Given is "digital information" in the form of a finite sequence {c\ : A G A}, 
where the data are taken from a finite alphabet, usually just c\ G { — 1,1} or 
c A G{±l±z}. 

2. D/A conversion: To transmit these data, they are converted to an analog 
signal of the form 

/ = $^c M 7r(^)#. 

AtGA 

for a suitable pulse g. Usually the lattice is taken to be a rectangular lattice A = 
ctL x /5Z. Clearly, the coefficients must be uniquely determined by /. Therefore, 
and in contrast to the previous discussions, the common assumption in wireless 
communication is that Q(g, A) is a Riesz sequence (for a closed subspace) in L 2 ' 
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3. Transmission of the signal f: The analog signal / is transmitted by a sender. 
This is a physical process subject to the laws of physics, in particular the wave 
equation. 

4. Distortion of f: During the transmission the signal is distorted by various 
effects. The most common effects are time lags due to reflection at obstacles and 
the Doppler effect due to relative motion between transmitter and receiver. Thus 
the distortion can be modelled by a superposition of time-frequency shifts, and the 
received signal is of the form 

(22) f(t)= [ a( V , U )M v T_ u f(t)dudr ] . 

Here the weighting factor a models the physical details of the transmission, such 
as the reflectivities of obstacles and the Doppler effect. It is well known that the 
distortion (1221) is precisely the pseudodifferential operator with symbol a, thus 

(23) / = a(x, D)f = ]T c^(x, D){ir(jJL)g) . 

Under the natural assumption of a maximum Doppler shift u Q and a maximum time 
lag r , the weight function a possesses a compact support in [0, r ] x [— z/ , u ]. If a 
were a bounded function, then the distortion operator a(x, D) would be a Hilbert- 
Schmidt operator and could not be invertible. Therefore Strohmer [20] proposed 
the generalized Sjostrand class M™'^ with an exponential weight v (C) = e a '^' as a 
suitable symbol class to model the distortion of time-varying systems in wireless 
communications. This class contains the distortion free channel corresponding to 
the identity operator (with symbol a = f ) and time- invariant channels correspond- 
ing to convolution operators. 

5. A/D conversion: At the receiver we decode the original digital information 
c\ from the distorted signal by taking correlations with time-frequency shifts as 
follows: 

y x = (/>(%) = ^c M fa(g,D)(7r(^)g),7r(A)g) < = (M(<t)c) a . 
AteA v 

The output vector y = (yx)x&A is completely determined by the received signal /. 
With definition ([9]) the input-output relation can be written as the infinite system 
of equations 

(24) y = M(a)c 

In wireless communications the matrix M(a) is called the channel matrix. 

For (!24j) to be well posed, the Gabor system Q(g, A) is assumed to be a Riesz 
basis for its span H 9i a = span£(g, A), and the distortion a(x, D) is assumed to be 
an invertible operator on L 2 (M. d ). Then the channel matrix M(a) is (boundedly) 
invertible on £ 2 (A), and the solution to (|24"|) is well defined. 

6. Equalization: To solve for the original information c we need to solve the 
system of equations f[2~4"|) and find the solution c = M(a) _1 y. 

At this point occurs amazing assumption that is taken for granted by engineers 
in this field: the channel matrix M(a) is assumed to be a diagonal matrix. 
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With this assumption, the distortion / is simply 

f = a(x,D)f = ^((7(x,D)7r(A)p,7r(A) 5 )(/,7r(A)(/)7r(A)p 
AeA 

= 52M(<rU(f,*(\)gM\)g. 

AeA 

In other words, wireless communications works with the implicit assumption that 
the distortion operator is a Gabor multiplier with respect to a Riesz sequence 
Q(g,A). The main motivation for this assumption seems to be its convenience and 
simplicity. With this assumption the solution to (124j) is simply 

c A = M(a)-\y x . 

One can show that the channel matrix with respect to a Gabor frame can never be a 
diagonal matrix. However, if a G M^ ' 1 with exponential weight v, then the channel 
matrix is almost diagonal with exponential decay. Theorem [8] then guarantees that 
the distortion operator a(x, D) is approximated extremely well by a finite number 
of shifted Gabor multipliers, or equivalently, M(a) is approximated extremely well 
by a banded matrix with few bands. In this case one may hope to improve the 
accuracy of inversion of (1241 by using a banded approximation of M(a) instead of 
the main diagonal only. In other words, we use an approximation of a(x, D) by a 
finite sum of shifted Gabor multipliers. 

This idea has been tested and implemented in [15,21]. Combined with other 
tricks, the approximation of the distortion operator by a sum of Gabor multipliers 
contributed significantly to a low complexity equalizer for time-varying systems. 
This equalizer performs best in strongly time-varying environments. 
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